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The modified Thomas algorithm is employed to obtain 
the correction vector AX in the widely used Naphthali- 
Sandholm method of solving separation processes problems. An 
efficient algorithm for obtaining the correction aX is 
proposed in which the sparsity of the submatrices of the 
Jacobian has been expoited in the matrix multiplication and 
inversion. The operation count for the proposed algorithm and 
the standard matrix multiplication and inversion has been 
presented using two bench-mark problems. It has been shown 
that the use of the proposed algorithm results in considerable 
saving in the CPU time with increase in the number of compo- 
nents and stages. The proposed algorithm can be used for the 
extensions and variants of Naphthali -Sand holm methods. 
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CHAPTER 1 
INTRODUCTION 

The classical techniques of separation of multicomponent 
mixtures like distillation, absorption, extraction, etc. ere 
widely used in chemical industries. The design or simulation 
of these (stage-wise) separation processes involves the solution 
of the material and energy balance equations and equilibrium 
relations for each stage, and requires an enormous amount of 
computational effort. Before the widespread use of computers, 
short-out methods were employed for the design though these 
are generally inadequate for systems other than the ones for 
which the equilibrium relations are linear. But with accessa- 
bility to powerful computers, the rigorous methods got impetus 
and several methods have been proposed. 

The rigorous methods of design or simulation involve two 
major steps; namely formulation of basic equations and their 
numerical methods of solution. The methods of formulation of 
the basic equations, in turn, can be classified as component- 
wise grouping of variables and stage-wise grouping of variables. 
The earlier methods of solution like the B-P method, the sum- 
rate method and the relaxation technique can be visualised 
as the direct substitution methods. Recently, these methods of 
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solution have given way to the more efficient Uewton-Raphson 
method and its variants. 

Under, the component -wise grouping of variables, the 
basic equations can he obtained as 

C ± Xi = - I i for i * 1, 2 ,...C (l) 

for the steady state conditions; where is the coefficient 
matrix involving vapor and liquid flow rates ara.d the k -values, 
\ could be either liquid or vapor mole fractions or component 
flow rates, are the feed component flow rates and 1 the 
components. In the relaxation technique (also known as Raise- 
Transient method) the basic equations are cast as 

= k+1 k+1 k , . 

CL X = X. i = 1,2... .C (2) 

i x 

where k is the iteration number. 

More generally, the basic equations may be represented as 

A . X as S’. I as 1,2... . .C (3 ) 

ill f 

The square matrix has the tridiagonal structure and can he 
found using the well known Thomas algorithm. 

Equation (3) together with the enthalpy balance 
equations around the stages oan he solved by the direct 
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substitution method for any assumed liquid and vapor flow rates 
and temperatures. The B-P method, the sum-rate method and the 
modified relaxation technique fall under this category* 


Since the year 1965, the direct substitution methods 
gave way to the Newton-Raphson method(or its modified versions) 
of solution. In these methods, the material balance discrepancy 
functions , 



and the enthalpy discrepancy functions around each stage are 


expanded in the Taylor series. Assuming the second and higher 
order terms to be negligible, the equations are rearranged to 


obtain 

J « b (5) 

where J is the Jacobian matrix of 2Nx2N. 

Several convergence schemes and the methods of obtaining the 
correction vector have been proposed. These have been discussed 
by Holland (2) • 


Naphthali (4-) and later Naphthali and Sandholm (5) 
showed that the convergence characteristics are better if the 
stage-wise grouping of variables together with the Newton- 
Raphson method are employed. Per each of the stages the 
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variables are the component liquid flow rates, temperature 

and the component vapor flow rates (i.e. l^,T,v^). The 
component and the enthalpy discrepancy functions are formulated 
to yield 

7 00 - o 

Applying the Newton-Raphson technique we get 


J 4 X « -7 

The Jacobian is of the order (2C+1 ) U x(2C+1 )N and it has 
the block tridiagona.1 form. Several methods proposed for 
solving the separation problems deal with the techniques of 
obtaining the correction vector 7 Stadtherjj) , 7, . 


It has been generally accepted that the Naphtha! i- 
Sandholm method has better convergence characteristics. However, 
it is likely to diverge if the guessed component liquid and 
vapor flow rates and temperatures are far from the correct 
values. To over come this problem, Ketchum proposed the 
fusion of the relaxation technique and the Naphthali-Sandholm 


method and demonstrated its suitablity for even interlinked 
columns. Later Hofeling and Seader r ij have demonstrated 
as to how the modified Thomas algorithm can be used for the 


interlinked columns. 


Stadther 



has compared the 
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convergence characteristics of the modified Thomas algorithm 
and the other methods of solution using the sparse matrix 
technique. He finds the Thomas algorithm as good as the other 
sparse matrix techniques. 

From the literature, it appears that the sparsity of the 
submatrices (diagonal and its adjacent submatrix) has not been 
exploited. The objective of the present 'work is to propose a 
more efficient method of solution taking advantage of the 
sparsity of the submatrices. 

In Chapter 2, the method of formulation of the "basic 
equations and the method of solution is presented. The results 
and discussion are presented in Chapter 3. Next conclusions 
is presented. 
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CHAPTER 2 

BASIC EQUATIO NS A HD METHOD OP SOLUTION 


In this Chapter, the strategy of formulating the hasie 
equations involving material and enthalpy balances, equilibrium 
and tray efficiency relations of multistage multicomponent 
separation processes, and the method of obtaining solution is 
presented* 


A model of a general stage, along with the notation 
employed herein, is depicted in figure 1 . The various discrepancy 
functions oan be formulated as follows: 

Por Material balance: 


M, = I . A T’ f, ~(1+S>. -(l+s)l. 

1 1 XL X p XX "* I X | Xl+ I X f XL H X f XX ^ X f XX 


(i) 


for l^ri^C & 1 ^ n <. H 


where v, and 1, are the flow rates of component i in 
i,n i f n 

liquid and vapor leaving the nth stage, and and s^l^ are 

the vapor and liquid side streams drawn on the nth stage. 


Por Enthalpy balance: 
0 

* 

E 


! n * Z'h.n-lYn-l + 2 T i. n+1 Hi * n+1 + y, 

*•' - <3 

- J> +S n>% ( n H 1(n - 2,(1%) li,n h l,n + \ 


n ♦ * ♦ 


for 1 f-i £ C & 1 X n X N 


( 2 ) 
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Equilibrium r 


elations together with the Kurphree efficiency 


°i,n 


^n k i.n h 


-t*£ + (1- \) 


i.n+1 


for 1 ^ i , 3 - 0 


& 1 4. n ^ N 


where 'l n = *^,**1 " ’ " 

Ihns, there are N(2C+l) set of nonlinear equations and they hare, 
to he solved to obtain the U(2C+1 ) unknown variables, m- 

1 . y . a nd 1 • 

i f n* i,n n 

The set of equations may he compactly written as 

f = F ( X ) - 0 (4) 


^"iitt. X , n+ 1_ 
■X- „ x, -J. . 


(3a) 


•where 


J 2 ^n 4 ”** ] ’ 

M . 0. # 0 o \ 

In **** C,n’ 1,11* 2,n n \ 


v 4 

II 

[ M 1 ,n ,M 2,n 

M 

•* ** C,n 

"V’ 

| f T . . „ , 

1 ~] 

JL — 

(JH * A 2 

_ > 

1 “i 


.In f 

ijH 

’ 1 2 ,n‘**‘* 

C,n 1, 


1 ( 8 ) 

nj 


It may he pointed out that the ordering of the variables is 

-from the one pronosed by Maphthali and 


slightly different from the one pr 
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Sandholm 1971 . The reasons for the order chosen here is 

explained later. 


Employing the Newton-Raphson technique, from Equation (4) 


■we can obtain 


^ X = 



where J is the Jacobian matrix. 

The Jacobian has the block tridiagonal' structure and can be 



where 

= ° n 

* = -r=- 

11 ^ Vi 



•A 




and the rest of the elements are null matrices. 
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In expand form, the sub matrix A^ is 


A = 
n 


i^.1 '5T 2 .n-l"“^ 1 0.»-1 ^ 1 ,n-r”Vc,r -1 ^n-1 


7>M. 


1 ,n 






*^2, n-1 




1 ,n - 


^J^n ,n , 




1 ,a 


C,n-1 ^ T 1 >**-1 ^ * 1 ** 1 ' 


> M C.n - * \ M C.n ^ ...^n 

^ 1 1,n-1 ^ 1 2*n-1 ^• L cX“1 ^ V 1,n-1 ^ v C,n-1 ^ n-1 



-1 

^0. , 

1 in-1 

>1 

,n-1 ^°1 .n-1 

^°1,n-1 


2^1 ,n- 

# • * * 

-1 

« o 

^ x 2,n-1 

« • 0 * « • 

• a 

,n-1 ^ v 1,n-1 

i v C,n-1 

• 

• 

^n-1 

# # * • 

?°C.n- 

» * 

zL 

^°C,n-1 

^°C.n 

-1 ^°C,n-1 

7o n . 

C,n-1 

?0 q^-i 

2 X 1,n- 

-1 

^ )1 2,n-1 

"c> 1 C,n 

-1 ^ ,n-1 

2"c,n~1 

^ T n-1 

d\ 


° n 

?*» - 

3E 
a n 

* B n - 


'd lj \ 

n-1 

fTT” 

° x 2,n-1 

^ 1 C,n-1 

*x rn ” ***" ♦ * ■* 

^1,n-1 

' > T 0,n-1 

^ T n-1 


ov ; 


0, 


V, 


0, 


h C °C C PL 


■where I, O f 0 are identity matrix, null matrix and null vector 
respectively. Subscript 0 indicates the dimension of the matrix 
or vector, h is the liquid enthalpy vector (h^ , i=1 ,C) and is 
heat capacity of I* ^ . 



fed li 
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In expanded form, the sub matrix EL is 



'JJM 

9 l.n 



"K.* 


9V 

3 M 2,n 


^ |1 1 ,n 
♦ * * # * 

’^2,n 

e * » « * 

'*♦### 

^ M C,r, 

^"c.n 

# « • * 9 

> M C.n 

Ti^ 

’^2,n 

...» *r-y- "* - t 

t> ± C,n 


^°1 ,n 


^l.n 

' ^ 1 2,n 
* * * # • 

3^C,n 

^°0jn 

'? r i,n 

• 0 * • « 

>C f n 

# • • * # 

^ °C.n 

••■•3 r c,n ’ 

^ E n 

3 E n 


^l ,n 


^ 1 G,n 


QM. 

^ 1 ,n 

^ M 1.n ^ M 1.n 

■3 T 1,n 

^ T C,n ^ T n 

^ T 1,n 

^ M 2,n *^ M 2,n 

$ T C,n 

# «■ * v * 


^C.n 

^"o.n ^“O.n 

^ ’ **» 

5°1,n 

2 v 1,n 

^°1,n ^°1.n 

2^’^ 

• • • • « 

• O • * • 

?°C,n 

■>1,n 

**•*• * * * 

3 °C.n ^°0.n 



^ E n 

3> 1,n 

^ E n- 3 E m 
2 T C, n ’ 3 T a 
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where X, X and X denote the filled matrix, vector and 
nonzero element respectively and the subscript C the orders 
of the matrix or veetor. 
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■^1 ,n+1 * 

>, 

_ii£ . 
^ 1 2,n+1 

^“l.n 1 M 1.n 

^ 1 C,n+1 f & Y 1 ,n+1-* 

2 M 1.n ^ M 1.n 

/ £ v C,n+l’ ^ T n+1 

7> M 2.n 

^2.n 

^ M 2.n F 2 .n 

^ M 2,n 

& v C,n+1 ^ T n+1 

2> 1 1,n+1 * 


^^0,11+1 ,n+1 

y * * * # 

f * «* * 


#*#*• • • ♦ 

W.n+l’ 

^2,n+1 

^ M C.n 

^^C,n+1 * ^ T 1 ,n+1 

^^C.n ^ M C,n 
^C.n+l’^n+l 

?°1,n 
,n+1 * 

* * * * ♦ 

3°1 n 

/■n 1 jJS # # # # # # 

,n+1 

3 °1_,n ^°1,n 

^°1 f n ^°1,n 

,n+1 ^ T n+1 

"»°0 A 

3o n 

C,n 

»°C.n > C .n 

^°C,n 3 °C,n 

5-^1 ,n+1 * 

^2,n+1 

^C.n+l’ 3T 1,n+l"' 

‘3'' r C,n+1 ’^ T n+1 



s>_\_ 5 !s_ 

^ v C,n+1 

^ 1 1 ,n+1 

^ 1 2,n+1 

^ 1 C,n+l’ 9 V 1 ,n+1 



o 11 
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unit "block 
submatrix 


null 

column 

rector 


null "block matrix 


Totally filled 

"block sub matrix 


**1 ,n+1 


C ,n+1 


and can be simplified to get 
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THOMAS ALGORITHM 

The correction vector 4 X can "be found by the Thomas 
algorithm and is given below 
Forward Substitution: 

For the stage 1 


Steps : 1 

P 1 -e 

— ( 





: 2 

Q 1 ^ 

--( 

h 

r 1 h 



For 

the 

stages 

n from 

2 'to 

(N-‘ 

: 3 

P <*- 

-( 

B 

- A P 

n- 1 

c . 


n 

n 

n n- 

“1 

n 

: 4 


— ( 

B 

n 

- A P 
n n- 

V' 1 

( F ■ 
n 


For the last stage N 


: 5 Qjj < ( B n - A n P n _ 1 ) 1 ( F N - ^ n Q n _ 1 ) 

Backward substitution 

: 6 2S X^j. Qjj 

For the stages n from N -1 to 1 

: 7 Al — -( Q - P X . ) 

^ n hi n n+1 1 

It can be seen that the algorithm involves several time taking 
matrix multiplications and inversions. Advantage of the sparcity 
and the stru.ctu.re of these matrices can be taken to minimise 
the machine operations, as described below. 
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In step 1 and 3, we encounter the matrix multiplications 

( B ) -1 0. and ( 1. -A.. P. )**<}.. The inverted matrices are 
* 0 _ I 3 _ »., v 

totally filled while C. and C. are highly sparse. The structure 

* t) 

of the resulting matrix is 







* 

-- rr 



Fill ed 


°C 1 0 ^0 


°c x o 5 o 

matrix 


0 0 x 0 5 C 

=r 

°c x c s c 



Ic Y c x - 


5 C x 


( B )‘V ( B - A } P..,) 0 d 


■where X, X and X denote the filled matrix, vector and non- 
zero element reepectiveljr, 0 and I are the null and identity 
matrices, and the subcript C the order of the matrix or vector. 


It can he seen that there is no need to compute the 
elements of the first C columns of the resulting matrix and 
we can save 0 (2C+1 operations. Further, the next C columns 
can he found as 


T ^ 

13,0+1 = b 3,i+c + * 

K-Otf 


'o,k %C+i 


for 


1 2C+1 

1 £ C 


(3) 


The elements of the last column are 

j 

P j , 2C+1 = * j , 2C+1 C 2C+1 , 2C+1 


(4) 



Thus the operation count of this matrix multiplication is 
(C 2 +l)(2C+l) instead of (2C+1 ) 5 required for standard matrix 
multiplication. 
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The matrix multiplication A P. .. (= c(.) carried out 

J J » J 

taking the advantage of spare! ty as given below . 




*c 


' ' t 

;°c x c 

*0 

% 

= 

O 11 
O 

0 il 
iO 

01 

o 

s 0 


°0 Jo 1 

*r>-l 


K n 


(5) 


Here, only the X^, X in the row need to be computed. 
The nonzero elements of are: 


^j.c+i - P 3,0+1 


for 1 3 4. 0 + 1 

i 

Uidc 


( 6 ) 




2C+1 ,C+i - 


^r' A 2C+1,k- 1,0+i + A 2C+1,20+i P 2C+1,0. 


for 0 + 1 (7) 

2 

Thus we need to perform only (0+1 ) scalar multiplication 
in stead of (2C+l)^. 
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The product A^ 

Vi 

si 

h-i 

. — 













*“c 


sa 






°o 

°c 

% 


*0 



J?0 


C HJ_ 


1 3? 

L 

f 

— — 


I 


3? 




■where 


P 


fi 


20+1 


3 

2a 


for 4 0 


20+1 > k 


\ + A 2C+1,2C+1 F 2C+1 


( 8 ) 


(9) 


only the last element need to he computed; thus only C+1 operations 

o 

are required instead of (20+1 ) operations for the standard 

matrix multiplication. In backward substitution, the multiplication 

P X .. involves (2C+l)(C+l) since the first 0 columns are 
n n+1 

zero elements. 

Inversion of ( B - A P . ) by partitioning. 

n n n-1 


Inversion by partitioning does not save the operations 

to be performed. But ftom the structure of the matrix 

( B - A P . ) some Saving in operations cah be realised as 
n n n-1 

(l+s n )l c X 0 

X C ^0 *0 

X 0 X 0 X 


shown below. 

(5 


n 


A P , ) 
n n-1 
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The inverse is given hy 




The inverse of can, in turn, he found hy partitioning 
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where 


- 


- 

- 


e 11 


+ 


e 22 d 21 ^ 


e i2 


•1 

a i2 e 22 



e 21 

= - e 22 

d 21 

toll )' 1 



e 22 

CM 

CM 

11 ^ 

11 

- d 

21 (a i1 )_1 3 

12> _1 

It may 

be : 

noted 

is 

a diagonal 

matrix and its inversion can 


Toe obtained by simply taking the reciprocal of each of the 
diagonal elements. Thus to find the inverse of a matrix of 
2C+1 , we need only the standard inversion of matrix of order C 
and some matrix multiplications. 



+) II 
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Operation Count 

(C 2 +1 )(2C+l) 

= 2C 3 + C 2 + 2G+1 

(C+1) 2 
= C 2 +2C+1 
(C+1 ) 


P X 2C 2 + 3C+1 
n n 

Inv 5C 3 + 14C 2 +6C 


3 

3 

n 


(N-1 )(2C 3 +C 2 +2C+1 
N ( C 2 + 2C +1 ) 

N ( 0 +• ) 

X (N~l)(C 2 + 2C+l) 
n 


Standard, 

operation 

( 20 + 1 ) 3 

= 8C 3 + 12C 2 + 6 C+1 
8C 3 + 1 2C 2 + 6C+1 


Saving 

6C 3 + 11C 2 + 4C 
6C 3 + 11C 2 + 4C+1 


(2C+1) 2 4C 2 + 3 C 

= 4C 2 + 4C+1 

4C 2 + 4C+1 2C 2 +C 

8C 3 + 1 20 2 + 6 C+1 3C 3 -2C 2 +1 

(N-1 )(8C 3 +12C 2 +6C+1 
N (8C 3 + 1 2C 2 + 6 C+1 ) 

N ( 4C 2 + 4C+1 ) 

(N-l)( 4C 2 + 4 C+1) 


Inv N(5C 3 + 14C 2 + 6C) 

Total - 3H(5C 3 + 15C 2 +9C+2) 
count ” o 

+ 2 (N-1 ) (2^+20^30+1 ) 


N (8C 3 + 12C 2 + 60+1 ) 

3N(l 6C 3 +28C 2 +1 6C+3 ) 3N(11C 3 + 1 3C 2 

+ 2 (N-1 ) (8C 3 +1 6C 2 +1 OC+2 ) + 70+1 ) 

+2 (N-1 ) (6C 3 + 
14C 2 + 7C+1 ) 
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CHAPTER 3 

RESULTS A HD DISCUSSION 

The algorithm proposed in the previous chapter has "been 
implemented in FORTRAN 1 0 on the DEC 10 system. First, the 
saving in CPU time for inventing the matrix (B - A C ) has 

XI Xl XI *""* I 

been examined* The efficiency of the proposed method has been 
tested by solving two absorption problems and the details are 
presented in this chapter. 

It is known that the inversion of matrix by partition 
does not result in saving the CPU time. But, the structure 
matrix in these separation process problems is such that we 
need to invert only C x C matrix to get the inversion of 
(2C+1 ) x (2C+i). This is so because the other submatrix, 
obtained on partitioning, whose inverse to be found to found 
is either identity matrix or a diagonal matrix with only a few 
elements which are other than one. Thus some CPU time saving 
can be achieved. As shown in the previous Chapter, the 
operation count for the inversion for a problem involving C 

f 3 2 

components is (5C + 14 c + 15c) compared the standard inversion 
(i.e. either by the Gaussian elemination or the Gauss -Jordon 

•X 

elimination) which is (2C+1 ) . 
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TABLE 2 

SAMPLE PROBLEM MO.I Taken from £ 3 ]} 

Absorption column having 20 plates and 4- component system;: 


Component 

Rich gas 

V N+1,i 
mole /hr 

Ah sorting 
liq. 

^0,1 mole/hr 

A 

75.0 

0.0 


B 

15.0 

0.0 

Absorbing liq. temp. = 125°C 

C 

10.0 

0.0 

Entering rich gas temp. =200 °C 

D 

0.0 

100.0 



Find the temperature and flow rates of the components of 
vapor and liquid stream leaving from the column. 



TABLE 3 


SAMPLE PROBLEM Uo .2 
STALE MEET OE THE PROBLEM 


C OMPONENT 

RICH 

GAS 

V N+1,i 

(mol/h) 

Lean 

Oil 

1 o.,i 

(mol/h) 

OTHER SPECIFICATIONS 

o 

o 

r\D 

0.4703 

0.0 

T q = 2.9°E, W 0*F, 

N 2 

0.1822 

0.0 

N = 8 , and 

CH 4 

88 .7000 

0.0 

P = 800 Ib/in^ abs 

C 2 H 6 

6.6747 

0.0 

initial temperature 

°3 H 8 

2.7786 

0.0015 

profile to he constant 

1 ° 4 H 10 

0.6375 

0.0006 

at T. = 25 °F 

3 

n C 4 H 10 

0.3655 

0.0013 

for all j (3 = 1,2...N). 

1 °5 H 12 

0.1158 

0.0067 

The initial vapor rate 

n 0 5 H 12 

0.0505 

0.0061 

profile is to he constant 

°6 H 14 

0.0146 

0.1495 

at V = 90.88(3=1,2. .. 8 ) 

D 

C 7 H 16 

0.0081 

0.5756 

and the liquid rates are 

C 8 H 18 

0.002 

1 .8214 

L. = 6.3095 (3 = 1,2. ..8) 

J 

°9 H 20 

0.000 

1 .6866 

and Lg = 15.42 

°10 H 22 

0.000 

2.0619 
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To get actual CPU time savings, numerical experiments 
■were carried out using the test matrix of the- type that are 
encountered in the separation process problems. The test 
matrices are tabulated in Appendix (D) . The CPU time required 
for the inversion by partition and by the Gauss -Jordon 
el emanation have been compared for different values of C and 
presented in Table 1 . 

The saving in CPU time is becoming increasingly 
significant as the C increases if the proposed method is 
employed. 

An absorption problem given by Naphthali-Sandholm 3 
and another absorption problem given by Holland 1 have been 
chosen for testing the efficiency of the proposed method. In 
the first problem, the total number components of the mixture 
are four and the number stages are twenty. In the other problem 
the number components are fourteen and stages are eight. The 
details of the problems are given in Table 2 and 3 and the 
equilibrium and enthalpy data are given in Appendices B & C. 

The elements of the Jocabian matrix (the partial derivatives) 
have been evaluated analytically. The analytical expressions 
for the derivatives are given in Appendix (A & E). 
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The two problems have been solved exploiting the 
sparsity of the submatrices and employing the 'standard ' matrix 
operations. Since no approximations were made while taking of 
the sparsity in computation, the number of iteration required 
identical in both the methods. The CPU time required to solve 
the two problems by these two methods alongwith other details 
are given in Table 4. 

Table 4 shows a substantial reduction in CPU time is 
achieved in the proposed method. The saving in CPU time becomes 
increasingly significant as the number of components increases. 
The proposed method can be employed even interlinked columns. 
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CHAPTER 4 
CONCLUSIONS 

An efficient algorithm for solving the separation 
processes problems by the -well known Naphtha! i-Sandholm method 
has been presented. In this algorithm the sparsity the sub- 
matrices of the Jacobian matrix is exploited in the matrix 
multiplications and in the inversion of the matrices 

(B - Ac ). It has been shown that the operation count for 
n n n-i 

the proposed algorithm is (5C^+ 14C^+ 6C) compared to (2C+1 

with the standard matrix operations. By solving two 'bench-mark' 

problems, it has been shown that the saving in the CPU time 

becomes increasingly significant as the number components 

involved becomes large. In the proposed algorithm, the saving 

in CPU time in computing the correction vector 4 X is effected, 
of 

but the number^ iteration required is the same as with the use 
of standard matrix operations. 

The algorithm can be extended for solving interlinked 
column employing the method suggested by Hocling and Seader. 

The computer code need to be tested to determine its effectiveness 
for distillation of nonidea.1 mixtures, extraction and adsorption 
problems. 
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N OMENCLATURE 

Subraa.tr ices of Jacobian Matrix at the jth row 

Submatrices of Jacobian Matrix at the jth row 

Submatrices of Jacobian Matrix at the jth row 

Total number of components involved 
Overall molar liquid flow rate from jth stage 

Overall molar Vapor flow rate from jth stage 
Molar liquid flow rate of component i from jth stage 
Molar vaprr flow rate of component i from jth stage 
Distribution coefficient of ith component at jth stage 
Vector of variables 
Jacobian matrix 

Molar feed rate of the component i into the jth stage 
Residual vector 

Residual vector for the 3th stage 

P matrix in the Thomas algorithm 

External heat input into the jth stage 

Molar specific liquid phase enthalpy of component i 

Molar specific vapor phase enthalpy of component i 

Molar specific feed enthalpies of component i 


Total number of stages 

Dfess balance discrepancy of component i for jth stage 
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s 


j 


Equilibrium relation discrepancy function for 
component i at jth stage 

Enthalpy balance discrepancy function for the jth stage 
Er action of the vapor stream withdrawn from the jth 
stage 

Eraction of the liquid stream withdrawn from the jth 
stage 


Greek letters 


difference in variable 

murphree efficiency of jth stage 

. Kronecker delta function 

10 

7) partial derivative 

activity coefficient of component i 
«_/V Thermodynamic parameter 

Zr P matrix of the Thomas Algorithm 

A. 

~n q matrix in the Thomas Algorithm 
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TABLE 4 

CPU time for two different problems 
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APPENDIX A 

ELEMENTS 0 E THE SUB MATRICES 
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APPENDIX B 
SAMPLE PROBLEM No. 1 


mm 

K values 


Material 

Temperature 



1 00°E . 

200°E 

A 

500.0 

550.0 

B 

1.50 

1 .8 

0 

0.90 

1 .00 

D 

1 .0 x 10" 6 

1,5 x 10-' 


Molar liquid enthalpies, 

10 5 cal/mole 

A 

0.01 

0.013 

B 

0.30 

0.33 

C 

0.40 

0.44 

D 

1 .50 

1 .90 


Molar vapor enthalpies, 

10 cal/mole 

A 

o 

o 

« 

v- 

1.002 

B 

1.80 

1 .82 

C 

2.00 

2.03 

D 

5.75 

5.95 
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Eor nonideal system only the following elements need to 
evaluate from different expressions, following expressions 
are evaluated using UNIQUAK model. Because of the lack of 
time anti for model independent programme, this feature was 

not incorporated. 
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Appendix E continued . 
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4 = exp(~ g-jjjp ) where a__ are UNIQUAK parameter. 
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Appendix E continued. 




MULTICOMPONENT 


s5sss====Bs=sasssssssgssssassps^asM-as*sssss*asssssssss.ssB 

DEVELOPED AND PROGRAMMED' BY L. PISSK 

nuwuff **®****®f*®**«> B S!*ssassSsa* ; s*assBSe*sa:KSJ s *iB 
EXTENDED NAPHTHALI SAND HOLM METHOD* 

SPARSITY AND SYMMETRY EXPLOITATION 

listing or symbols are as follow^ ^ — *— s= — 

s,*,s»s3**r***Bw*»»ajs#**sjSR*»3*»*a!ff«ff*jB afl 5* j i« SE B3 5 * S s. s ,ssass S gs S; igRe. 

AL(J)SLIQ, RATE OF THE J. TH STAGE’ 

vujjvIp rate of the. j th stags 

svcj.n ; yap, rate' of component i at the j th stage 

|o fiS S HEAT I ?ftg8| T lT 3 rHi 3 J P TH E s|ASSi^ r ***' 3 T “ STAG€ 

SSt J) .FRACTION OF VAP, STREAM TAKEN FROM THE J TH STAGE 

!?f{|i!«Sp| I E g\?f 1 p?4riri5l.^rg: !, i? ft II E J « STW£ 

Kill gijjftfg g f cB3?33iUr I Si 

P * IrlPfe v : rn* ?I?jR l E lla.I? 1 T P?H s ?£fif 0F IHE - TAGC 

I?nr»T * f Tft TU Dl*S- F r»J H SftbI82SSl Cll . r ^* !,, ®^- I S® 0 FIRST PLATE! 

SYoilJ? |Ap8p FLOW RATE OF COMp.ONEif r pSNTERING I liTO P N A Ta PLATE 

5§DE:« ! coiS§igfe R X^lgg MIE un ™ 1)113 "* * IH «-»« ■ 

.9L component i next op. stage of the stage 

, i Kl Xh,j§ I kI B 5il?S!l COMST ^ r valos; for component i at any stage 
n: number or stages 

cs TOTAL NUMBER OF COMPONENTS INVOLVED FOR SEPARATION 

S£W«iIiES E S-Slii.8RJ D !!S 3 SSS T , i: .ftF. I ?S : . J TH SfhQE - 

frc JJJFEEp TEMPERATURE OF THE J TH! STAGE 

0 : ENTHALPY balnce normalisation factor 

S = .SSSSSSSCSSSSSS9SSXSSSSSSSS3SS^ ; S,S i S.SsS3SSBS£SSSSSSS=SBSS^BSSS 

SUBROUTINE ENL ; COMPUTES LIQUID* PHASE! ENTHALPY OF ALU 

COMPONENTS FOR A! GIVEN STAGE 
SUBROUTINE ENV * COMPUTES VAPOR PHASE! ENTHALPY OF ALL 

COMPONENTS AT AS GIVEN STAGE. „ „ * « 

DIST ! COMPUTES THE) EQUILIBRIUM CONSTANT FOR ALL* 

COMPONENTS AT’ GIVEN STAGE 

DDIST S JOM|uIes THE: DSp.VA|IVE OF EQUILIBRIUM. 

CONSTANT K W.R.TV TEMPERATURE FOR AS S l AGE 

• • V ■ DEWL J COMPUTES THE DERIVATIVE OF THE LIQUID: / 

PHASE ENTHALPY OF ALL COMPONENTS FOR A 
PARTICULAR STAGES 

»'# V. DENV J COMPUTES THE DERIVATIVE. OF THE VAPOR .PHAgP 

ENTHALPY OF ALL. COMPONENTS FOR A PARTICULAR 
STAGE 

BI l COMPUTES THE INVERSION OF B SOBMAIRIX 

^ bi i : Compute inversion of bi matrix 

MEN VECTORS 

RIEGER C.CT 

COMMON /ApEM / SLO»SVB!» TB.T0, VB/AR£A2/5Sf SS5, IF ,ETA, SQ 


i * '. * ■* 
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DIST 

DDIST 


DENL 


........ 



itflMlilti 



SEC 14, 4) 
TF(8), 


Hft (45 u 

4ErA(a)#AK(i4).r 
5 VO ( 20) * SC (14,4) 

OPEN OJMIT*2l, FILER'S, DAT') ,, * s *' 

;S;;?;::*m;SS5? S w?^?;: s 2 sss r ssssaasaaa * 3a3ss * 3s3 »=s~= 3 ««s a3 

§d £S5 VECTOR’ wo* DMA SPECIFICATION 
s ! 8 f s ? st f! sSsasssssss * s! * s * ses *i l!s;s ** s SsW B s s *«sssassss 


sa.satssaSSB#. 
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READ 
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READ 
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READ 
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READ 
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11 :! 
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j{l 



r J*; » ?! 
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Jisl. 
Jfcl,N 


N) 


C),J*i, N) 


mmm 


READC41,* 

READC41,* 

READ(41 ,* 

mm:i 

RKE(4!:*5i«ii!:5Kssi:ij;isi:»« 

END OF INPUT SPECIFICATION 
sstxasscsassanstssssisssszs 

AEF*0. 

701 

FORMAT? 4X, 'INITIAL 1 GUESS OF FLO* RATES AND TEMPERATURE! ARE: ') 
H Ft X T fP C 4[ 4k % 2 H 1 1 ' **" ' 

forhw(/;JL*vapor phase: component flo* rates are*) 

WHIHl»|jlj».M) 

FoJSAT(Jjc?^MPONEMr NO. MQX, 'STAGE NOS, ARE') 

m i* I 

FQRMATC26x jt M\12X r , 2M2X, , 3M2X,*4 , ,12X, , 5' f 12X,'6 # ,l45X, 

1 * 7 %i 2 X, * 8 f ) 

NR1T|C44:224|) 

F0RWATC4X f 116C *-*)) 

fpH»r|/>j4X®‘i,I8. — - 

8811111:!* 
mimUV'i 


WRITEC44, 

NRITEC44, 



COMP.' FLOW RATES' ARE' ) 


.'INITIAL VAPOR AND Limo RA 


fOR^JS//, 'INITIAL VAPOR AMD LIQUID R 
1 ARE # // ) 

WRITE(44 f 4467) 

WRlfl 14:7979$ 


2HID RATE. AMD TEMPERATURE GUESS 


WRITEC44.6100)' 

JgRMAT(//^OX f »CO«POSei»T *,SX,'ro?, Wff,fEEO\ tOX, 'BOTTO«i VAP,F E £: 
F3RMAT?|6x|io| 4JC. 18 ♦»!')) 

fl»;p lHki 1 

NRITECS, 8 ) 

WRlTiCi,H) 

NRITEC 5# 15 $ 

F0RMAT(/,4x# 'DEVELOPED AMD 1 PROGRAMMED' BY U. FIENK’) 
**a#**s*a**a*8ss»asai!sssssa8sssMj*?s*«sass5sssssas*as*!S?ssSB 
I COMPUTATION OF AUGMENTED ERROR FUNCTION AT THE GUESS VECTOR 
| sssasaaaeaaaaasajsaaassaaaaaa.s.saaaaaaaaaxssassaasassasa.saKSa*. 

# # # CONTRIBUTION DUE. TO MASS BALANCE ###' 

sasassssssssassssistasssssssssssiiassjsssrssssssssssssss-asss*. 
EVA^UATlON^g| L |Oa) & V(Ntl r I) FOR 5 TOTAL CONDENSER 

3s:ss3s3ss33s:3s:ss3sss»sssssaa3a8aas:as8asssss 

8 nD*qP CONDENSER & REBOILER PART 

IF ENTERING LIQ.. AT PLATE 1 IS NOT AT TCI) ADJUST THIS TL' 

accordingly 

DO*"l Jel , N-l 

:kp i, . 

igcjjNf.i) Goto 7$ 

iPlD^I 

?om||nu1 

FlCN?I)I-h£cN-l,nfSVB(I)+SFC«,lO -Pt»SVCN,I)-p2*SLCN,I)) 

CONTINUE 

C2SI^!^4iiiIL^^J^J25,ppSSSiJiWij. s 

I END OF EVALUATION OF MASS BALANCE CONTRIBUTION 

i Sts3lSS3*3S3*BB*«SSIS3SS38M38SS-S«2-*S-***** 

I START OF EVALUATION OF THE! CONTRIBUTION DUE TO 

# aaaa«assaaaaraaaaaaaaa*a».*aaaaaaaaeaaaa8aaa8aa,. 
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£35 3 » 
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E*|+St-Cfi-i;h*SHH(I)fSVBia.)*HH(I3fSFC»,X)*HF(I3-P2»SIiC«,i3*SH(I} 

h(N,CT)»-(E+5Q(N2) 

CM# C.T3 


( Ui 


*) 


jWWRfliStc,. , 

WRfTEC5,22)N,jle 
FPF ® ht*W 

« TO COMPUTE THE LIQUID & Y&P3R FLOW RATE FOR &Lfi< STAGES 

# ssa8ssassss*i?sa*«*8a«*.ssss?ss3**saaspsssss5#j?«38«s, 

aaaxaKeaaasasasasasaaaaaaasasaa.ssaaasBSBsasasBsssBsassssBSjs.sets. 

« ?g??g5|PJ°3r D ^L I lri® !,WRWM ' c ° NOmON co»siotRi«G 

I a*aa*#a*aa»»»aaaa*aa»ala«a«a«aa*a»*'B»aa*aaaa*aaaP«aae, 



TO 4 JslfW-l 
rr*,T( J) 

AtiL DIST(TT,*K) 
>0 4 J*l,C 
i a# *t* X 


J(5+1 {)/ifjtH ) ) fll!C(I)$SHJ,I)/:!VWJ3 * $V(J#U/VCJ)+cl * s - AcJi)) * 

^lGr{ri(j;c+i5»Fi(j,:^i) 

^NTINUE 
TT*T(M) 

CALL DlSTjTT,AK) 

m5j|+I)*^CETA ( N)*AK; c I)4St ( N,lO/AUN)-5VCN,I)/V£N) + c l.-mCN)) 

iEF*AEF^FlcN f C+I)*Fl (M,CffIO 
CONTINUE 

33s»:c3z?s:3ss3asssss3:s;ssss::3£:s?i;s::::s3ss3:ss:s:s: s sss. 

#. END of error FUNCTION EVALUATION 

s.S333:sssss2$sssz:$3ss::s:s:3:t,ES353£7:s:£ssr=::::s3s:ss33s3. 

IF(FFF-EPS) 5,5,6 ' 

a*,*»#*aK*»**»«s*»#!»x«st»»sxjsjBSfc»s,3«,a ( saa?a?sKBSs:s=a=SB=psss3ssss« t 

FIND THE DIRECTION OF THE N IH PLATE; 

«« rT*T(L 

pi«4.+§§ 

P2»I ,*SS 
TALL PIS 

:all ene 
:all en‘ 

‘all od 

>Q 1$ I 



rPA) 


5§VS^‘ 

mm 


)®PA(I) 

¥ ,T,PA) * 

im 

ALL"0|NV(TT,PA) 

Oo J* 1 a?,,. 

If ljC)*PACI) 

‘ ,/AL( 1 ) 

■0, 

VV**j/V?l)-2 

p;||;g|n,i^8(ifi) 

88 a! iSi ijc 


s tm, «# 

!JJ Sl 

tilk 





g£8S=»sSSSS*;c*SiS c ^' a ^ 1>w ^ wpas 

Du 4§ I«1 »C 

ISJ'SIt?. Pi 

8v I rC JIl* "Pi 

pD OF ^COMPUTATION OF Bi 

k staiEt’ 

HH f T ) *0 

)+SV(2,I)*PA(I) 


4, 3X) ) ) 



|T|NUE 

;ia+c,n+c)*(i*Etxci))»(vc2}"3V(2! f i.))*P3 

,QNTI | • #r i 

j it jh » q n t i*h.zn 



HETACi ) ) *5VC2f I) *P3< 


fe P ?a!feS '« raK ' 

mil l\t 

§Ui 

w . - ,cn, !*.!,:!) 


/;(4X,9CEU t 4,3X))) 



w 

S?Ir| ( 0#' 3 ?N^ C AND^WWO^STAGeS 


ERROR* SECTORS / C4X,9tE13. 6 »3X))) 


iUUtihl 

P2at,4§s|(5) 



ENt(TT,H) 
r,SHH) 

feiH5 Ha) 








TION 0 




C 5 , 30i ) I (Cl (I , Ji ) , JI«d , CD ' 1*4 cm 


mmw 


W tis? 

iimuw 

CALL EMV(TtIh) 

5 AM, DDI5TCTT, PA) 

do &o M,e_ % 

?Si5tHuk 

CALiti DKNVCf T, PA) 
;m.L de«l(ti,p») 

feiMiSUi, 

M N 

/hfo*2 


•)*0 

;* ' 

1 iy 

i*»P 
[)«ipi* 
'*B{CT, 
__ _u(J»I)* 

I0NTIMUE 
)0 85 1*1, C 
>3 85_n*Cc 


FCX),I«1,C.T) 



i»8(C:tI,CI)*P3*SKJi,I.) 

)io 


•sv 

I) 









#pltTT*P2 
*0 


C?LL*S§SM*j 

SHI'S ft* 

«s*ni 5 iii 


ioi I# JI) , JI*4 ,CD , 1*4 # w® ) 

c 6? i ; ji ) , ji*i # c n, i*4 f cm 

i,ji),ji*4fCf),r*4,cr) 



,3X))) 



DO ifO Ipl.GT 

"iff 1 f t > 

fc&iriSufc \ 
k« m t5» Gnro 70 

THESE A ELjUfiMTS ARE UEHEMTS OF Z\ SUBMATRIX 

$$M SShhcd'o 

m 9 D kT' pk) 

?5 mtinue 

wr Bl 


s lUlsSM 58 is 1 1 :=*?* H » S »* , .» > 

DO 5? HJJ i{* 

y ih.xfWum,,!) 

*Rir£(5 r 'ioi)((cia,jn,Ji»l.cr),iM,oi) 

?h!*r|ti;if-ra) 

MrmU? 4x H # }aL 6 es a£?i M 11-1 ' * 

p f i ?fcp iiwLi-wc« 

IVklikBr 

c* S* JS°* 

CONTINUE 

55 ii 

|;jH» 

p ' x “* c ** 

M,C* 




tjftjli •sv(j; i 5 -oNrH 3 ;!fC) 

ln55r oF,F1(J ' CT) 

PoImATC/', 'ItIrATIOM HUMBER* M3) 
AEF=Q. 

03 75 J»1,N 

[ftklsg* 

B8 13| 

: 0 NriNUE C ^ ( i * ii>ib0 * 



IBITE(44,9) . 

iT(/,:§UCESSFOt, CONVERGENT') 


F3RM|\T( 


:u 

jRAMON NUMBER* ',130 

IAL VALUES OF FLO* RATES' AND TEMPERATURE* AXE J * ) 


r 1 0 X , ' S I AGE: MOS. ARE') 


PHASE. COMPONENT* FUON RATES ARE * ) 

||b., ... 

fn§An 2 {>lill%l 2 K,’ 2 ' , 16 *.' 

17 ,i 2 Kf 'ft '5 


f f m \* \ \ 

4,41815 
//) 

M * W * * 



, M) ,I*1-,C0 
$ fO 

FLOW RATES ARE'' ) 


I , ( Shi Jif I) # SWi ,H), 1*4 •#■ CD 



VAPOR AND LIQUID RATE AND TEMPERATURE ARE 


4X i 


FORMA|( 4 S: A 42 r.*)) 

lit !i 

IFAIL =1 

STOP 

eSd 

SUBROUTINE MTMUL2(A,B) 

§S 3 f ,, ?I??c| fc, '= T> ' 1 

B * 1 lS, 68 T a jm,ci 

9 UrJ}«D(I,J) 


B(CT,CT) rD(2t>29) 


8&i*»5S- 

ETURN 


Sy|R5UTINE_MT«ULlCAf 8) 

,B(CT),D(29) 
i*B(II)+DCI) 


ypnyuu-wfe ^ a rou mm % it 

MTEGER C # CT . 

imensiSg a{ct,st),i 



iUBRDUTINE ENL(T,EL) 

i^S|i#^C),SCCU, 4 , 

CONTINUE 

RETURN 

END 

SUBROUTINE DENLCTi EL) 



ELCC) r SC (14,4) 

>0 : i ) *si 1 5 , 2 ) + 2 . *sc ( u 3 ) * * * sac t, 4 ) *r *2 

ITINUE 


LIQUID 



return 

^ROUTINE ENV(T,EV) 
* N 


COMMON Cj 
jvjMQN /A 
[PENSION 


C),SE(14,4) 


I * 1 ) +$E( I , 2) *T+SE'(I,# 3 ) #T , *2!fSE{ I*4)*T*3 


t* t y 1 1 c" 

kM«nnll ' ^ 

RETURN 

fm 5 

SUBROUTINE DENV(T,EV) 

INTEGER c,ct 

COMMON C.CTrN 

DIMENSION EV?C) ,SE(14,4) 
n n i IssIjiC 

EV(I)sSE(I , 2)+SE(I # 3)*2.*T>3,*SE;(i:#4)*T*2 

CONTINUE 

return 

IuErOUIINE 5 UBMAT(B,Cm 

m^7c C t C r\ 

''1£NSI56J{CT,CT),C1.CCT,CX) 

1 1*1 rCT 

BU,J)-C1(I,J0 

m r 



SA(H,4) 

n+SA(i: # 2)*T+SACl:, 30*r*2*SAa, 4)*T*3)*3 

U0RQOTINE DDI8TCT, AK1) 

i NfESER CiCT 

ESmon /ilt/iA 

Pension akcc),sa(14, 4) 

il 

CONTINUE 
RETURN 

INF°GBf?Bf RAt * 8 MATRIX NIT* SCARCITY EXPLOITATION 

mmm 1 9 

Dimension bJct, CT),BSC i4r 14) # asu i4,H) #PA( t4) f PB ci4) 

pA(h**l>i(i f i) 


#T*( 

*r 



J)a6§lCl,J)+BU+C,IJ:)*BS(IJ,, J.) 


OMTINUC ^ 

„0 2 1 * 1 r C 
DO 2 J*1#C 

gMitilS*" nwI ' ,+a 
B8 1 S=i:S 

P l m?‘ 

HLL FMI(BSl) 

0 5 t*l 



yUNJTANUl- _ 

DO 6 r*i,c 

DO 6 J»Lg 

SSVfjii 

"Ms 

DO 7 J*UC 

Ifcx.nUici) 

*\hhh , 

! # j)+bs(i,J) 

9 

vri| 4 wy& 

8B*BB+S(6f,I)*PA(I)tB(C!rrI+C)* ! PB(I) 

i(:T?C|)»l./(BCCTrCn*B80 


Cl) 






*ajs3ssasssaaE»a*iaaES3Ss:ssss=s*3!asas^*,ssjs«,*,s5 1 sssssssaB 

SUBROUTINE Emu) 

INTEGER C»CT 
COMMON C,CT e N 


i C) , TT ( 14, 14) , UI( H, 14) , T(14f 14) ,ETA(14, 14) 


B »L/A(K,K) 



DO ' J 1 r ^ 

7 1»ui(I,J)tETA(IrK)*UI.(lC r vJ) 

JUE 

1 = 1 ,C 

_ , J»l ,c 

IP |sB° 8 

go 9 j*k{£ 

grj|iJi)ii?i,u)*Aaji,j)+ETMi,J?) 

§9 }i ljU S 

|a yEiifiw) 

IFCK.EO.l) GOTO 5 


88 !! 3 = 1:8 

DO It XJftl f C 

I! hfc 


§nn*v 


RETURN 

END 

INVERSION 01 MATRIX 

DURING INVERSION SPARCin IS EXPWirtO. AND COMPUTATIONAL^* 

FASTER THAf^AS USUAL INVERSION 


SUBROUTINE BII(B) 

coP|S r c c £t T n 

DIMENSIOA^B {CT# CT) f SB ( 1 4 1 1 4) # P Al(14) # PS ( 14) 

PA(I)»1. >6(1.1) 

PStf )*pJU I+C) 

CDNIINUE ' 

DO \ 1 *} , f 


H » ir 

I f Jj *0 # 

CONTINUE 
DO 4 t*l,C 

i(lTC<*U-B(I+Cf J)*PAS(J) 

B( I . J J«0. 

CONffMUE 

Igcit ,a)+-aci» u+so*B(iJ*ct, J» 

shHsj.j 

contJnJe 0 * 

|8 l| SijS 



BJTSSt.c „ 

8 CI+C# J)«SB(I, J) 

pMffNU? 

DO 10 1*1 eC„ 

PAcl)=PMl 5 *BCI, h 

pbuwbuS+bu+c, 
bbsbb+b( 6 t , 


fimiftij^iftggissjaas^cT) 



n*PA(i)+B(::r,i+cmaci.) 

,al./(BCT,Cr)-BB) 


gi? II fUl c c 

CONTINUE 


J+C;) 


14 1*1 

5r, x+O* 
Mtinuc 



i.) 



3K j ) 


'll) 


JUTINE ACMUHArU, 

5pIi|iSfi^cT#cir)#ci<Cf#Cf) 


i mf- 

>o 1 5*i ii 

i»£l(l,X)+ACCT,J)*C*<'M+CD 

:+i)*ci(i,i)+ACCTrCf)*ci(crri:*ca 

<tJ « ¥ fniir f#I) * clcI,cr) 

!i(Cr,CT)»AA+ACCTrCT)*C:l(CI#CT) 

ftAaO. 

RErtJ&N 

E*D 


nmnh 

SlMRNSIofi AfCT,CT),FCCT) 
po i t*i,c 

F(:ff!!AA+F(CT)*A(CT,:nP) 

REruftN 

lu^RpUTINF BCMULt 8 , C 1 ) 

jN 

00 2! 

00 2 J*l,C 
DO 2 K*1#CT 

gl(f . J)*ct CI # J)fBCl,K!)*:i(K,J+:i) 

8“ ! 

if:-" 

^sitPiE I#t)4,BCI,K!) * ca(K ' CT) 

OOJ t = l,CT 

| fi s - 

return 

B e f% 


l)^B(I»K:)*C3l(K»CT) 


»iilS I c E ci FMUt,tcl,r) 


B fesvAk „ 

NSl§fici(CT,Cr),FCCf) 

1 * 1 ,Cf 

ji)|6T(i,i)tci(i,KtC9«rcctici) 


DIMENSION 

00 t 1st, 



h ’ 


tr [;v 

-rrino 

tip 

PRO 

1 ** 0 1 

,27 6- 

03 

.06 

4 *'*m U ? 

,626- 

0 3 

.93 

7 7 p m» j i i 

i ;i?> 

02 

, 9 l 






, Mr-0 2 


,«9 


. !’■>' -*'2 . 1 

, 4E-OJ .99 

* ■■■■■■ <■- * 3. 

.55{>02 . 49 

9 2, M <? 

» ■ ' •«’ # ,f 

, 7E - ' . 3 

■* 5 P 1 

.926-0? , yp 

♦ r; r* m ix r ) i 

m ■ ' ■ w •», 

.11 Fl-Q 1 ,88 

. 1 , t 

. 15E-01 ,8« 

. ‘ ,1 

. 10R-UJ .87 


*5 


' r - iV , 
• A j 

* 45^2 
’ 


^ r t ^- * * 9 \ * • f? ' 

W-Pi .g|E*0j ,37 
96E-01 . 17E-01. .36 
?4K»0| .4SR-01 .35 


. - .91R-0* ,56£»01 . „ 

. ’ '>"«»? .60E-OJ .34 

. M«fi-02 .85E-01 .32E-01 .83 

. 7,>f *>02 ,8|r+02 ,§6E+02 ,R8E+Q2 .89E+02 .9E+02 .91E+02 
.'AEtO§ * n J E+02 .92E+02 ,92F+0? ,926 + 02 ,976+02 ,036+02 
, ARf 02 ,946+92 ,956+02 ,966+02 .976+92 .936+02 

: j f JS2, * * •. ' E !•?.''• . J 1 * f ? J . * 1 » 6 t ?.? «»»*♦•* 


.1 ”•« 


nr+''3 

t I !'+ •■< 3 


’TY s 1 ’■ r -T“ * * * * otv 5 .ijrt ij , l . 
11K + 0 5 . 1 I 6 + O 3 .116+0 3 . 1.1F+0 3 
llE + i* i .11R+03 ,116+03 ,116 + 03 


:! ? 5 
I /* ’ 


l 5 7 

1 *> t. , 1. 9 4 

» ; _ r 


? r >7, 1 6 0 , 
197, 


16 3* 

200 * 


166, 170* 173* 176, 


c; r i , ' 

5* 

4 • 

7*4, 

7' * 


6 * 7 ’ 

a 


t, (■■■ fi i r 

• • s ' 

. " J * J • * 

1 * f 5 * l 9 « 

t ^ "V • ^ * 

! 6 • 1 9y . 


6,9 5 

6 

6 


■ 4 


' 0 

r* 

' n 


i \ 


. ( 

». i 

J 7 

0, 

), 


0, 

* 

0 . 

0 

• 

. 0 

1 . 

1 * 

1 . 

1 

* 

: 1 

0,. 

' o « 

«. 

0 

• 

0 

M 

■A 

■ + ’» 

n 


■ ; 0 

* 

» 

*■ * 

• 


V « 

O' * 

■T’; 

a 


0 

■0 :■ 


f \ * ■ 

* 

" 0 


■. 5 

«» 

' * 

• 

' .0 

# 

(j 

0 * 

0, 

f? ■ 

* 

■f) 

* 

.. A 


1 

o 

6; 


0, 


o* o 


0 * 0 
, 0 
. I 
0 


0* 0* 0, 


0* o, 
0 . o , 

0* 0, 

■ r t> . 


0. °f 

o .5 


* u *V * V * 

0. 0. 5, 0 
o. o. o, .o 


. 0, 200 

0, 0, 

i - i. 

0, 9, 

0. 0, 

0, o, 

5 ,5 



MftiKHor 


RATE or 

putai. r» 


S5!!f5S^^ 0 l§uJ l ,?ggiS8 s IS?^S?iiSWigS c f "mncwrowM 

f552*?SS55*™**f *•*«**•»»»** **»«**, >***«. pxaamgmuiiBMBmtisttstBJ****- 

PROGRAMMED BY L,< FIENK 

a ,»«* »«» a««sa»*a **»«**»* wo aa*|t»sE^«,s l 3s*s***E«5*jeaaES:sscss^<?<!c 
NAPHTHALI SANDHOLM METHOD 

*»*» bJE«J*b*i» *«**»**» wasass* a* 8*B*a»J8Sff«a*^s>s*=eas*a#BSjKsas^« l 
LISTING or SYMBOLS ARE AS FOLLOWS' 

*nxmuM9*M*m*mat*9*Km*xm*i**nn*Mxi!tm#mi*8,p*nmm9**X9wa!*s('*Bae;pFB‘ 

v( JO i vap Rate or the j Th stage- 

syij,i)i Vap, rate or component' i: at the j th stage 

lkH^Uai I ?A 0 RrL 5 LS^ p 5 S e gJ«l : ^ T vm J TH stage 
so j heat Input at the j th stage) 

§§S?n F *IM2f«4 F J A ri* s I8fa2 TAicp from the j. th stage 

Ilf? J? I ^RgACTf QN Or L|Q w irftiAM rgiREH FROM THE J TH. STAGE! 

STAC JJiMURpHPEE^ErrlCllMCY^gr THE J! th: STAGE 

SH(l5f LI0 PHASE, ENTHALPT OF Cfinio»g||' I AT ANY STAGE 

HI |5| VAP, PHASE ENTHALPT OF COMPONENT I AT J TH STAGE! 

HFfljj ENTHALPY OF FEED AT ANY STAGEi FOR COMPONENT’ i: „ „ 

lI'fHjglS^^iJLCilfPpNflir I jIlLfJTTOM STAGE OF THE: STAGE 
TB I TEMPERATURE OF THE VAPOR ENTERING. N TH PLATE 
TO i TEMPERATURE' OF THE LIOUID ENTERING- INTO FIRST PLATE. 

SLO(I) I LID, RATE Of COMPONENT i: ENTERING AT FIRST PLATE. w 
sy5(l)» VAPOR FLON RATE OF COMPONENT' I. ENTERING- INTO M* TH PLftlE 
VS f TOTAL VAPOR RATE ENERINS. INTO THE N TH PLATE: 

UNDER CONSIDERATION 

Mliii™ i c * ,Mp3Nc ' Mr 1 M£xr up! sTftGE ° f m sTftsEi 

AK(f?| KOI JlBRIUM CONSTANT VALUE* FOR COMPONENT I AT ANY STAGE 

Nt NUMBER OF STAGES 


?S ( 5 

TO * 


r fSiH mr 


AK(I)I EOILIBRIUM CONSTANT VALUE* FOR COMPONENT I AT ANY STAGE 

S! NUMBER Sf Sf AGES , 


SUBROUTINE CHI. SSDHRUTES MfiUIU' PH»S|, ENTHRLPT Of UU * 


SUBROUTINE ENV JOOMPUTES VRPOR' PHRSEi ENTHRL 

COMPONENTS ITJtJim- STJP 


* 0IST *cSMPOMeirs r fcf : SHrSttK CDNS " NI F0R “ LLI 

DEHO , co|| pyipWfc iMfeM * 

<»"* . =p|l«>iplI f ^5^g e 5?s r ?gR^ P ??R??f§^R 

: th 

i «* m m ns SB me m m as M m M fi 2 S 2 S M S £ JS 21 *£ 2 & 22 P # 2 22 * 2 22 *■ 2 2 2 2 2 2 22' 2- 2 22 2*8 ’ 




SSW8 8 HUti/su 
R?3Sfe,.£M 

DIMENSION B(29,29 


» T0,VB/ARS] 
29),PB(14 


, TF , ETA , SQ 
C2C8, 29,29) 


'Mr 



SE(14,4) 

TFC8), 


5V0(20),SC(14,4) 
a,****Fi»*a*B**»**iB**i*«i«»*a»*****jB#*j«jta#t#*»}8|B#*t*ai^BjBs****i«!l*esis, 
OPEN (UNIT*21,FILEa , S*0M') 

INPUT TERMS MAINLX G4IESS VECTOR AND DATA SPECIFICATION ' 
uMMM »** M as ini S 3 ft m M *ft® ftftftft z ttM sjsftzsssijBxssss^ss; mm. mft. 

#*}N f CT*C T1 T . , ... a , 

.*! h¥ fU " 


m\tu* 


m 

READ 
REAP 
READ 
READ 

m 

READ 

READMli* 
READ (41 , * 


READ f 41 # * 
REAOl 41 ' * 
READC 41 ,* 
READ (41 *-* 



- j»i , n 


C),J*1,M) 


J \ \QR,% l t <J J ijmX , l&l 0 wQ 

Il*SSSS«Kllt8a«SS9XS«8XXS8SX:CStM5S!ISWS«SSSS«S;B8S?. 

SND : OF input SPECIFICATION 

SSItZSSSISsXBXtaCSlSBStlSlB 

ww mm upF mP mpf p®' mm Mppi mm imp 'iw ; mm mm mm mm mm fm* -mm mm mm -"mp imp * 

' ~F*0. 

0 701 

TIALt GUESS OF FLON RATES AND TEMPERATURE} ARE? '5 
Hapor phase: COMPONENT* rum rates are:') 

^ ,, 

ponent NO. M OX, 'STAGE. nos.- ARE*) 
),10X,14('«*),/5 

12X,'2M2X,'3M*X,?4Si2X,'5',12X,'6',m, 





m*,co 


COMP*f FMJW RATES' ARE * ) 


tSL(Jl,I),JW#N)#lM4v CD 


VHP OR AND RUfe. AMD TEMPERATURE GUfiSs 

MIlliM.????! .. .... _ 


FEED ',1 OX, •Borrow vap.fee. 


iN6 " ■> 

mml 

FORMAT!?, 4X» 'DEVELOPED AND* PROGRAMMED* BSf 


mlmrm 

FORMAT! ?,4jL 'DEVELOPED AND* PROGRAMMED BSf L, FIENK ' ) 

• sjtaasasattaftsPaaaSsaaaaaaaaasaaaaaaxaPPPsssPSPsssssaaaa.spaaae 

I COMPUTATION OF AUGMENTED ERROR FUNCTION AT THE GUESS VECTOR 

# 5*B*atwa8B*ssaas3#*ss*ssss*s,?* l *«?a:MW»aP=sass?sss3s5f^=sse s - 

9 » # CONTRIBUTION DUE, TO MASS’ balance: # * # 

aasBaxacassasssssasataasaasaaaaapaaaaaaaasassaaspsxaasssspasaaa. 

EVALUATIgN B 0|^0(I) A V!NH,I) FORi TOTAL- CONDENSER 
a****!#»a*jEaB**a*aaa»s»eaaesa*a*aaaaaaasa»*s*»*ss 

ffcy?K8 E Z?^, IrTdm tt'Sw AT m> adjust ihts » 

.ACCORDINGLY 


^KRiSr* 1 ”'' M 

ssthmii 

i : k -Ui r. 78 


SRntinOe 
GOTO 1 


,I)tSF!Ji,i:) HP1*SV !J,I) -P2*SL( Of I) ) 
SS BALANCE*.', El 1.4) 


f | C I ) if { ?J S l| ** 1 ' 15 .P1*SVC 0* I) -P2*SL»( Jir I ) ) 

i^Li# f ?! f I j ?nf lisf i ? ‘ ^ ^ « - p * * sv c “ * 1 > ** ,2 * si - <** ' 1 ^ 

illlm l l ' ll 

PirE(44,18)AEF 


END OF EVALUATION OF MASS' BALANCE, CONTRIBUTION 


STAR? OF EVALUATION OF THE! CONTRIBUTION DUE TO. 
E£l2SSt£t«M*a*£ia*aaaB*aaaaaaa*»l»ia»«P*»*»*S!»'«*. 


hums 


K#«M*J(KXSK«Sj*tt»asa*K,XJ 

E *0 . 

do 5 
fr*.T(j) 
pl*l. + SS(J) 

P2*'I,+SSS( J) 
rFF«TF(J) 

Rft 

$p:ir&s?& 26 

Fi n cl ) 1 ' 

GO TO £112 

g »3 

pfp|p|i|jSS??5«tS 


!tb,hh) 


jssexjc.; 


i#x«p»x«xxx*xxx*x«*sxa. 




I)*HF(I) 


ft rt T t 

' w f A 

r I Jnil 


1(1, CT) 

4X,'fcEF0F EWmiiPr 






f f S ^fSl f J i * ??} X > + S F c a , I , . «r t ! , 


fimi 4 


CJirCD 


S|oL ENt* 
CALL ENV 

e*- 5 . 

DO J.72 I 




) +S VS ( I ) *HH C I) * SF C N , I ) *HFCI)-P2*SLCN, 1 ) *S« C I) 


SOniMUg '7 

WRffi?J 4 5i?J N AEF r 

# TO COMPUTE THE 

# X*BX»****X**'*»X.**»*****J | P*-*«S t «?*^*** ,s|! ' , ‘* a * ,!=l 5 SS **C*“f c f ,: 

J g?PfgfiS?l o gr°XSc, ir IrlSiI^ BRI0M CMDITI0M CMW ^ ,fi 

#. »*EiSiSS«««»i»#»wM*****»*««*^*«**»»***»*«**!»?«- 

mm*mBmm9*mnmn*9mmu*mmmm*%n*m*Miz***9mii9m*****9******* l **? t ** 9 *‘ 


jB*a»xx- 






DO 4 J*1 »N»1 
CA*| C oisX(TT#AK) 

CALL IIsTCTf i AK) 


Ua+i l) 7 vfj{Kr RK]Cl)fSLCJ,I)/ ^' ( j 3 * Sv ^#I 5 /V(J.) + (l-EaA(J 0 )* 
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